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Abstract 

This paper develops the theoretical foundations for the ability of a control field to cooperate with 
noise in the manipulation of quantum dynamics. The noise enters as run-to-run variations in the 
control amplitudes, phases and frequencies with the observation being an ensemble average over 
many runs as is commonly done in the laboratory Weak field perturbation theory is developed 
to show that noise in the amplitude and frequency components of the control field can enhance 
the process of population transfer in a multilevel ladder system. The analytical results in this 
paper support the point that under suitable conditions an optimal field can cooperate with noise 
to improve the control outcome. 
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I. INTRODUCTION 

Control of quantum processes is actively being pursued theoretically [1, 2] and 
experimentally [3, 4]. In practice, control field noise and environmental interactions in- 
evitably are involved. The present paper considers the influence of field noise upon the 
controlled dynamics with the noise described by shot-to-shot pulse variations, as in typical 
signal averaged experiments. Recent studies considered several aspects of the influence of 
laser noise[5-ll], and this work aims to further explore the issue. The interaction between 
the noise and field driven dynamics is generally a highly nonlinear process. The impact of 
the noise can either be constructive[12-15] or destructive [10] in the manipulation of quan- 
tum dynamics as well as either improve[16, 17] or reduce[18] the convergence rate of control 
search efforts. These disparate behaviors make it difficult to precisely identify the role of 
noise under various circumstances, but the many successful experiments at least support 
the point that some noise can be tolerated [19-24]. Operating under closed-loop [25] in the 
laboratory will naturally deal with noise as best as possible. A theoretical analysis[9] on the 
impact of field noise upon optimal control indicated that an inherent degree of robustness 
can be anticipated by virtue of the controlled observable expectation values being bilin- 
ear in the evolution operator and its adjoint. Some simulations of closed-loop experiments 
also show that robust control is possible and that properly designed control fields can fight 
against noise [8, 11, 16, 17]. 
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Recent numerical simulations [11] of closed-loop control in a model system showed that 
under suitable circumstance the control field could cooperate with the presence of noise 
to more efiiciently reach the target state. The cooperation was remarkable when seeking a 
modest yield (e.g., ~ 10%) in the target state. In this case the noise or the deterministic field 
acting alone would each produced a small yield, but the two acting together cooperatively 
produced a much larger yield. The latter numerical simulations did not reveal the underlying 
physical origin of the cooperative efi^ect, and the present paper will analyze the controlled 
dynamics of a multistate system in several limiting cases to explain the prior findings. The 
cooperative behavior also has foundation in analogous stochastic resonance phenomena[26- 
28] and fluctuation control[29, 30]. 

Section II presents a general control model of population transfer in multilevel systems. 
The noise is modeled by run-to-run variations in the control amplitudes, phases and frequen- 
cies with the observation being an ensemble average over many runs as is commonly done 
in the laboratory. The goal of the control is to maximally populate a highly excited state. 
In Sec. Ill we develop a weak field perturbation theory to provide an analytical solution to 
the outcome of the controlled dynamics. We obtain the noise-average yield from applying 
the control field in Sec. IV. It is shown that variation of the field phases from pulse-to-pulse 
plays no role in the dynamics but strong cooperation between the deterministic portion of 
the field and noise, in both the amplitude and frequency, is possible. Finally, we draw some 
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conclusions and discuss general multistate systems in Sec. V. 



II. THE MODEL SYSTEM 

The effects of field noise on controlled quantum dynamics will be explored in the context 
of population transfer in multilevel systems characterized by the Hamiltonian H, 

H = Ho-i^E{t), (la) 
Ho^^Enln) {n\, (lb) 

n 

where \n) is an eigenstate of Hq with the associated energy in the absence of radiation, 
and II is the dipole operator, fj, — Ylm n' l'^) ('^'l- '^^^ control field E{t) has the form 
which may be implemented in the laboratory [31] , 

M 

E{t) = 2s{t)^Ai COS {ujit + 9i), (2) 
1=1 

where {cc;;} are the frequencies of the radiation, and s{t) is the pulse envelope function. The 
controls are the amplitudes {^4;} and phases {9i}. 

Noise in the laboratory could take on various forms and arise from a number of 
sources[10]. In keeping with laboratory practice, the achieved control will be measured as 
an ensemble average over the outcome of many noise contaminated control fields, and the 
noise is modeled as shot-to-shot uncertainties in the amplitudes Ai, phases 9i and frequencies 
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ui inEq. (2) [8, 11]: 

A; = A° + X,, (3a) 

ei = e'l + yu (3b) 

uji = + zi, (3c) 

with {xi) — (yi) — (zi) — 0. The shot-to-shot field noises are captured in terms of zero- 
mean uncertainties x — {xi}, y — {yi}, and z — {zi}. For simphcity we assume that the 
different noise components are independent with distributions p\^^ {xi), pf^ (yi) and (zi), 
respectively. In practice the field noise in E {t) may have complex structures and origins. 
Besides the uncertainties in the control field amplitudes {Ai} and phases {Oi}, a potential 
source of additional uncertainties is in the frequencies {uji} due to laser frequency jitter from 
variations of the refractive index[32], or other sources [33]. The fiexible treatment of random 
variations in {Ai}, {9i} and {cui} is meant to represent all of these various possibilities. The 
net outcome of the control experiments is the average, 

O [E it) ' 7] = f p!^^ (xi) dxi j'\ pf^ (yi) dyi pf^ {zi) dz)j O [E{t, x, y, z)] , (4) 

where 0[E{t, x, y, z)] is the control yield 

O [E {t, X, y,z)]^\ (^/IV' m, X, y, z),T])\\ (5) 

produced by the field E {t, x, y, z) in Eq. (2) using the amplitudes, phases and frequencies 
in Eq. (3). The target state is |^/), and \'4> [E{t, x, y, z),T]) is the state of the field-driven 



system at the final time T, which is a functional over time of E{t, x,y, z), t < T. In what 
follows we assume E (t) ^ for t — > ±00. 

The objective function to be minimized with respect to and {Of} in the presence of 
noise has the form 

J= |0[E(t),7] -OTp + aFo, (6a) 

i^o = EK)'> (6b) 
I 

where Or is the target value, and Fq is the fluence of the control field whose contribution is 
weighted by the constant, a > 0. 

III. WEAK FIELD PERTURBATION THEORY 

To illustrate the principle of how the deterministic portion of the control field can coop- 
erate with the noise, we consider the excitation along a ladder (or chain) of nondegenerate 
transitions and energy levels with each linked only to its nearest neighbors. One could 
analogously think of the system as a nonlinear oscillator [34-36] or a spin with 5" > 1 and 
nonequidistant energy levels. The transition elements are taken to have the form 

(7) 

The + 1 level system consists of an initially occupied ground state |0)ati— >— 00, A^ — 1 
intermediate states \n), n — 1, 2, ■ ■ ■ , N — 1, and a final target state \N). The states are 
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coupled with an external laser pulse having the nominal form of Eq. (2). The wave function 
is expanded in the form 

N 

m-J2^n{t)\n)e-''-'. (8) 

n=0 

The initial condition at i — > — oo specifies that Cq — 1 and Cn — ior < n < N. The goal 
is to maximize \CN{t)f with \ijjf) — \N) for i — > oo, when the field is zero. The result of 
perturbation theory for Cn to the lowest order in the control field E (t) is 

Civ = ^^^'"M dtME{tM)e'^^'^ dt^_iE(t^_i)e^---i*--i . . . / dhE{h)e'^''\ (9) 

f.—]^ J —oo J —oo J —oo 

where 

U^n^ ^n- ^n-l.n^l,- ■ ■ (10) 

are the transition frequencies. Utilizing the Fourier transform of the field 

1 f°° 

E{t)^—j /(n)e-"*dQ, (11a) 

/oo 
E (t) e'^'dt (lib) 
-oo 

in Eq. (9) produces 

N „ 



^ k=l 

/oo ^ />iiv _ l't2 

■oo J —oo J —oo 



:i2) 



In the above equation, the definite integrals over the frequencies Jl^ run from — oo to oo. To 
calculate the time integral in Eq. (12), we introduce a small imaginary part — > +0 in the 







Nf{^l)---f{^N) 


exp 


—I 


Y^^=l {^q - ^q) 




-1 
1 





frequencies, 

Qfe^Qfe + i£fe,A; = l,--- ,A^-1. (13) 

It is seen from Eq. (9) that such a change will not affect the result, since E{t) — > 
for t — oo; we assume that E{t) decays at least exponentially for t — > ±oo, which is 
physically reasonable. Evaluating the integrals in Eq. (12) with respect to ti, t2, • • • , ^at-i 
yields 

TV 

/■oo _ exp -12^^^ y^Lq-UJq)lN 

■^-^ [-^T^,=l{^p-^p)\ 
where all Q^, k — 1, • • • , A'" — 1, carry a hidden small positive imaginary part as Eq. (13). 

Integrating Eq. (14) again with respect to ijv and Jljv produces Cat in the frequency repre- 
sentation. 

An important property arising from Eq. (15) is that for the state N to be populated by the 

pulse E{t), the sum of the transition frequencies should be equal to en — sq. There is an 

important difference with the work of Larsen and Bloembergen[34], in which the condition 

£n — ^0 — Nuj leads to multiphoton Rabi oscillations, not to an actual transition, as they 

have a stationary periodic field, not a radiation pulse. 
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We will consider the case where M, the number of the components in the pulse Eq. (2), 
is equal to N, the number of the transitions in the multilevel ladder system, so 



N 



E{t) = s (t) J2 ^/e^'^'e-^'"'* + c.c, (16a) 
1=1 

N 

f{uj) = [Aie-'''S {uj - uji) + Aie^'^S {u; + uji)] , (16b) 



1=1 



where 



/oo 
s (t) e'^^'dt, (17) 
■oo 

and each component is only resonant with the corresponding system transition, 

\uJk-ojk\ < \uJk-ojjjtk\ , = 1, • • • ,N (18) 

with the function S (a;) assumed to be smooth, \dS/duj\ <^ \S\/uJk. The latter condition 
means that the typical duration of the pulse significantly exceeds the transition periods 
2t^ /ujk, k — 1, - ■ ■ ,N. Then the nearly resonant terms in Eq. (15) are kept, 

J 1 _i l^p=i\Sh-^p) 



N 



= C'^J]/XfcAfee-^^^ (19b) 
fe=i 

where 

t^i j=i Ep=i i^p - ^p) 

If each component of the pulse is exactly resonant with a transition of the system, 

ujk = Ok: (21) 
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then a simple form for Cjv, 



C„=^, (22) 



may be attained with 

/oo 
s (t) dt (23) 
-oo 

being the effective pulse duration. It's easier to prove Eq. (22) in the time domain than in 
the frequency domain. Considering the control field in Eq. (16a), whose each component 
is resonant with a particular transition, the neglect of all nonresonant terms in the Eq. (9) 
yields 

^ poo l>tff l>t2 

CN^i'^YlpikAke-''^' dtNs{tN) o?ijv-is(i;v-i) ■ ■ • / dt,s{h) (24a) 

I 1 J —oo J —oo J —oo 



k=l 

N 



I 

k=l 

■N N ^ 
I T 



N 



r s(f)df 

i"" n i^'^'^^'"'" ' ^ m ^^^^^ 



YlptkAke-''", (24c) 



k=l 

from which Eq. (22) follows. 

If the pulse is not exactly resonant 

Sk^UJk-Oky^ 0, (25) 

then Cn is not so simple, except for the case where Sk — S is independent of k, i.e., the 
detuning is the same for all frequencies. In this case Eq. (22) still applies, but Eq. (23) 
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should be modified to become 



oo 

-iSt, 



S{-6)= / s{t)e-'''dt. (26) 

J —oo 

Introducing the following change of variable in Eq. (20), 

z^ = n^-Uj, (27) 

yields 

where is the cumulant detuning, 

k 

Afc = ,iV. (29) 

p=i 

Inserting delta functions in Eq. (28) produces 



/ iyV-1 ■ / /.oo /.oo /.oo \ N-1 _/^^-;' 

-T^^ivTiyn ^4/ c^r, 5(- ^ - A^.) n ' ■ 

(30) 

From Eq. (13) it follows that there is a small positive imaginary part in each variable in 
Eq. (27), so the integration over z'l^, k = 1, ■ ■ ■ , N — 1 yields 

■N ^-^ / /.oo r°° \ 

= 7^7^ n / E ^« - n G(^,)e^^^(^-^^^+^^). (31) 

(27r) fc^i V^o J-oo / j.^i 

The same result can also be obtained directly from Eq. (12) by changing from ti, • • • , In-i 

to Ti = ti, T2 ^ h — ti, ■ ■ ■ ,tn — tN-i (or ijfc = Ti H h Tfc, /c = 1, • • • ,N) and from 
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Qj to Zj given by Eq. (27). As illustrations of the general formulation above, Gaussian and 
rectangular pulses are considered in the following treatment. 



A. Gaussian pulse 



If the pulse envelope of the control field is Gaussian, 



s (t) = exp , 



S {uj) — rexp 



9 / ' 



(32a) 
(32b) 



with (7 being the spectral width of the pulse. 



a — 



2x/i 



(33) 



then substituting Eq. (32b) into Eq. (31) leads to 



C 



■N N 
t T 



N 



{2n 



k=l 



(27r 



oo poo 

dTk 
-'-oo 

oo /-oo 



dzk ] exp 



A 



N-l N-l 



N 



7V-1 



a 



j=l 3=1 \p=l 



7V-1 



N-l 



k,j=l 



(34) 



where the elements of matrix A and vector h are 



Akj — ^kj + 1; 

N-l 



bj = I 



E-.-2- 



^N 



(35a) 
(35b) 



p=j 
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Carrying out the Gaussian integrals in Eq. (34) with respect to Zk, k — 1, • • • , N — 1 yields 



C 



N 



N-l 



(27r) 



71 



N-l 



detA 



■N-l 

n 

k=l 



It's easy to verify that 



dTk exp 



N-l 



N-l 



^N 



k,j=l 



(36) 



{A ')^.^6kj--, 



detA^N. 



(37a) 
(37b) 



Performing some algebra produces the perturbative solution for the Gaussian pulse, 



C 



,N^N 



N 



■>N-l 



TT 2 



■N2 



2 N-l 



n 



exp 



exp 



Pn{t)\, (38) 



where 



Dk = kl^N - NAk, k^l,--- ,N-1, 



N-l 



PAr(f) = ^A;(7V-A;)T,^+ J] 2k{N-j)nTj. 

0<k<j<N 



k=l 



(39a) 
(39b) 



When Df; ^ and cr — > 0, the asymptotic behavior of transition probability [3 7], 



Cn\ oc a 



N-l 



exp 



Na^ 



(40) 



is obtained. The above equation shows that the yield exponentiaUy decays with the sum of 
the detuning of the individual transitions, An from Eq. (29). This behavior can be easily 
understood since, as pointed out previously, the sum of the field frequencies should be equal 
to the sum of the transition frequencies of the system, which is the energy conservation 
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condition. A simple calculation shows that, for all \Si\ ^ a, the most probable way of 



meeting the energy conservation constraint leads to 



C 



N 



(X ex.p{-A%/Na^). 



B. Rectangular pulse 



In the case that the pulse envelope of the control field is rectangular, 



1, 0<t<T, 
0, otherwise. 



lUJ 



the substitution of Eq. (41b) into Eq. (28) leads to 



C 



N 



T / ^iv-i n 



^-2 iziT _ I 



dzk 



{-2n) 

where In-i is an integral with respect to -Zat-i 



'-N-l 



-I 



exp 


-i {j2q=l + ^iv) 


T 




—i 







exp 


—i 


(j2q=l Zq + ^n) ^ 


^] -1 dZN^i 


—i 


T.q=l Zq + ^n] 





(41a) 



(41b) 



(42) 



(43) 

This equation can be evaluated by the residue theorem[37] . After checking the three poles 
in the lower half plane Im z^-i < 0, it can be shown that only the residue of the zero pole is 
needed because those of the other two poles do not contribute after integrating with respect 
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to zn-2-i so we have 



1n-i — 



-27riRes {z^-i = 0) 
2m exp 



EN-2 , A 
q=l + 



EN-2 , A 



(44) 



Integrating with respect to zn-2, zn-s,- ■ • ,Z2, similarly as above, reduces Cjv to an integral 



with respect to only one variable, 



^ _ {-lfi r e-^("i+^^)^ e'^'^-l dzi 
(-1)"^? /• e-^(^i+^^)^ -1 dzi 



M A ■ (45) 

This result finally produces a compact form for the perturbative solution with a rectangular 
pulse: 



Cn = / " , , , , c^^. (46) 



If all of the transitions are resonant; Sj, = 0, p = 1, - ■ ■ ,N, then 

(-1)"; fe-"'^ 



which is consistent with Eq. (22). If all A^, g = 1, • • • , A?", are different from each other, 
then from Eq. (46), 

TV AT 

CW=(-1)^E^"'^^^"^^^^ n (A.-^)"' (48) 
9=0 j=0UH) 
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with Ao = 0. In the case where all of the detunings are the same, we have = q5, then 

Cn = - if. (49) 

This expression agrees with Eq. (22) derived earher for a pulse of an arbitrary shape. 

An interesting consequence of Eq. (49) is that the transition probabihty is an oscillating 
function of T5, 



C 



N 



2 



(TV!)' 2 



In particular, for TS — 2mr, n — 0,1, • • • , 



C 



N 



becomes zero. Such an antiresonance is a 



result of destructive interference, which eliminates transitions to higher states. 

IV. COOPERATION BETWEEN A WEAK FIELD AND NOISE FOR A MULTI- 
STATE LADDER SYSTEM 

Prom Eq. (19b), it is evident that, in a weak field driven multilevel ladder system, the 
transition probability is independent of the phases of the individual pulse components, but 
is sensitive to their amplitudes and frequencies. This section considers the noise- averaged 
transition probability over the amplitudes and frequencies of a weak field, and shows that 
the population transfer can be enhanced under suitable conditions. 
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A. Noise in the amplitudes of the control field 



Prom Eq. (19b), it is easy to determine the dependence of the control yield on the 



amplitudes of the field pulses, 



N 



0[E{t)]^\Cr,\'^a%l[Al 



(51) 



1=1 



where 



N 



k=l 



(52) 



is independent of the amplitudes. If each amplitude Ai is contaminated with random uniform 
noise distributed over [—7/, 7^], as in Eq. (3a), then the noise- averaged outcome of the control 



process is 



N 



0[E{t)n]^a%l[{A^), 

1=1 

where (A^) is the mean square amplitude of the l-th radiation component. 



(53a) 



(Al) = J (A° + ^ifpl^^ i^i) dxi =Af + (x^) . 



(54) 



It is instructive to compare this expression with the outcome from the control field not 



having amplitude fluctuations. 



N 



0(°) [E it)] = a%l[ A\ 



02 



(55) 



1=1 
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The ratio 

1=0 

becomes appreciable for a large number of states N, even when the ratio (Af) / A° ^ is close 
to 1.0 for each individual transition. For example, if (Af) / A'^'^ — 1 + e with £ <S 1, then 
0/0(0) ^ 1 + 7V£. 

To find the optimal field amplitude, we set (^4^) = A^'^ + {xf) and minimize the functional 



N 



N 



+ aJ2Af (57) 



1=0 



1=1 

over ^4°. Assuming that are independent of the optimal value of ^4° satisfies 

Af + {xf)^ Constant, (58) 

independent of /. An important consequence is that, for given {{xf)}, the contribution from 
noise can beneficially serve to decrease the required amplitude A° of the optimal control 
field leading to a given yield, provided that the yield is small. Cooperation with noise can 
be extended to modest control yields beyond the perturbation approximation, as shown in 
numerical simulations[ll]. Although the presence of strong noise can considerably reduce 
the coherent nature of the dynamics, modest target yields can still be reached, including in a 
very efficient manner. However, when attempting to achieve high control yields, a different 
mechanism is generally operative involving competition between the deterministic portion 
of the control field and the noise[ll]. 
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B. Noise in the control frequency spectrum 



Here we consider a weak control field with a Gaussian envelope [Eq. (32)] and frequency 



noise having Gaussian distribution 



fk (h) 



1 



cxp 



4 - 4 

dk(T 



(59) 



where the random variables are the noise contaminated detunings 



h = 5k + Zk = ujl + Zk- iOk 



(60) 



with {zk} being uncertainties in the laser frequencies {cuk} in Eq. (3c). An additional 
source of the detuning noise can be from the transition frequencies {cuk} due to the Doppler 
or phonon-induced shift in the control medium. Recalling Eq. (38) yields the noise- averaged 



transition probability, 



q=l \J-oo / 



[Pn (r) + Pr, (f')] 



(61) 



Here, Lat is 



N 



q=l 



N-1 



^ ^N-^^Y.^rk-Tl)Dk 



k=l 



W n ( / d6Acxp - J26kBk,S, + J2hck-J2^: 



(62) 
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with the elements of the matrix B and the vector c being 



Bkj 



(63a) 
(63b) 



and the parameters T4 and V specified by 



3=k 



1 ^ 



(64a) 
(64b) 



Integrating Eq. (62) with respect to 5q yields 



nLi4VdetS 



■ exp 



■ 2 N 



with the elements oi B ^ (inverse matrix of B) and the determinant of B being 



(65) 



kj 



dl I Skj — 



2d) 



N{l + 2d?) / ' 



TV 



det B = ^, (P = — 2J 4- 



nL4 



AT 



fc=i 



(66a) 
(66b) 



The expression for the scaled transition probability is simplified in two limiting cases. 
The first case is that of nearly resonant driving, where \5k\ <^ adk. In this case = 
—ia~^ iy — Vk) and the term proportional to 81 in Eq. (65), can be neglected. However, 
it is easy to see that Ljv (t, r') < 1 because of the noise in the frequency spectrum, and 
Ljv rapidly decreases with the increasing noise intensity parameters d^. Therefore, the 
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transition probabilities (|Civ|^) with noise are smaller than without noise when the control 
field frequencies are reliably tuned to be resonant with the system transition frequencies. 

The role of the noise is reversed in the case of comparatively strong detuning, ^ adk- 
In this case we have fa 25k/(j'^d\, then Eqs. (61) and (65) produce 



CatI^) oc cxp 



(67) 



N(T^ (1 + 2(P) 

with Ajv = SfcLi ^k- In this region, it follows from Eq. (67) that increasing the noise 
intensity parameter leads to an exponentially strong increase in the transition probability 
This finding shows how noise can play a constructive role in a controlled quantum system. 
This result may be understood from the fact that for strong detuning, spectral noise can lead 
to some pulses actually being closer to resonance. For these pulses the transition probability 
is exponentially higher than for nonresonant pulses. As a result, the noise averaged transition 
probability is also strongly increased. 

Noise-induced enhancement of the transition probability occurs also for rectangular laser 
pulses. This is most easy to see when the detunings of all the frequency components in 
the pulse are the same and the scaled transition probability is given by Eq. (50). As noted 
earlier, the transition amplitude into the target state is completely eliminated if the detuning 
satisfies 5 — 2m: /T. Pulse-to- pulse variation of 5, or pulse-to- pulse variation of the duration 
T will suppress the antiresonance and lead to a nonzero transition probability even when 



C 



N 



for the average values 5 — 5,T — T. If the width of the distribution over 5 is small 
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compared to 5, but largely exceeds tt/T, then from Eq. (50) we have 



N 



(2N)\c-2N 



V. CONCLUSIONS 



This paper explores the dynamics of population transfer in a multistate ladder quantum 
system driven by noisy control pulses, with particular emphasis on identifying circumstances 
when cooperation between the field and noise may occur. The noise enters as run-to-run 
uncertainties in the control amplitudes, phases and frequencies with the observation being an 
ensemble average over many runs as is commonly done in the laboratory. If the rotating wave 
approximation is valid, the quantum dynamics in the weak field limit is greatly simplified 
and independent of the control phases. Furthermore, if the objective yield is modest, the 
control field can cooperate with amplitude noise to reduce the applied fiuence. Frequency 
noise in the control field is shown to be capable of enhancing the transition probability when 
the detuning is large. In the laboratory implementation of closed-loop control, the optimal 
field will be deduced to extract as much beneficial value as possible from the presence of 
noise. This paper presents a theoretical foundation showing that ample opportunity exists 
to take advantage of noise. 

The above conclusions are fully consistent with recent numerical simulations[ll]. Al- 
though the analytical treatment in this paper only applies to ladder-configuration systems. 
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the basic conclusions on the prospects for cooperating with noise should be applicable for 
optimal control of many multistate systems. This point is confirmed with a non-ladder 
system [11] where a high degree of cooperation was found between the noise and the de- 
terministic part of the field. Regardless of whether it is dynamically beneficial to fight or 
cooperate with noise, the optimal field will appropriately emphasize the dynamical pathways 
that correspondingly either work with or circumvent the influence of the noise[ll]. 
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